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Abstract: We propose a novel design of an all-dielectric optical antenna 
based on photonic-band-gap confinement. Specifically, we have engineered 
the photonic-crystal dipole mode to have broad spectral response (Q -70) 
and well-directed vertical-radiation by introducing a plane mirror below the 
cavity. Considerably large local electric-field intensity enhancement -4,500 
is expected from the proposed design for a normally incident planewave. 
Furthermore, an analytic model developed based on coupled-mode theory 
predicts that the electric-field intensity enhancement can easily be over 
100,000 by employing reasonably high-g (-10,000) resonators. 
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1. Introduction 

In analogy to the radio-frequency antenna, an optical antenna can function either as a 
transmitter (emitter) or a receiver in the optical frequency domain. For use as a receiver, 
incident electromagnetic wave's energy shall be efficiently captured by the antenna at a 
certain resonance condition, resulting in strongly-enhanced local-electromagnetic-fields. 
Recently, various metallic nano-antennas have been proposed and demonstrated [1-6] to 
obtain such strong local-fields (or hot-spots) for bio-chemical sensing applications [6,7] and 
for boosting a wide range of nonlinear optics phenomena [8,9]. However, the performance of 
these antennas is often limited by their low Q < -20, which is due to either huge radiation 
losses or absorption losses in the metallic structure itself [4]. Moreover, it has been known 
that resonance frequencies of plasmonic nanostructures do not have a strong dependence on 
scale [10], which makes it challenging to develop an array of diverse plasmonic resonances. 

In this paper, we propose to use a 2-D photonic-crystal (PhC) slab structure [1 1] as a novel 
platform for the efficient optical antenna with strong local field enhancement. Until now, 
various 2-D PhC cavities have been extensively researched for achieving ultrasmall mode- 
volume (V) and high-g enabled by the photonic-band-gap (PBG) in the in-plane directions 
[11-15]. Evidently, such non-metallic PhC cavities can provide much higher-g (easily over 
10,000) in comparison with their metallic counterparts. Furthermore, the resonances of PhC 
cavities are scalable with the lattice constant, allowing lithographic means to tune the 
resonances and to integrate them in a dense array within a small footprint. 

2. Coupled mode theory for photonic-crystal optical antennas 

Our scheme for a PhC optical antenna is shown in Fig. 1(a). A high-permittivity 
semiconductor PhC slab cavity is positioned nearby a plane mirror. By varying the gap 
distance, d, one can control the directionality of the far-field emission, which in turn will be 
used to optimize the antenna coupling efficiency [16]. An electromagnetic planewave with 

intensity I Q is being illuminated at a normal direction (# = 0), where the incident electric- 
field is polarized along the x-axis. To understand the underlying physics, we first consider a 
simplified model in the spirit of coupled-mode theory (CMT) [17,18] as described in Fig. 
1(b). The energy amplitude in the cavity denoted by a is being built up by the incident 
planewave while decaying into (infinitely many) multiple 'ports', which are imaginary two- 
way waveguides connecting the near-zone to the far-zone and vice versa. A symbol is 
assigned to the port lying in the direction of k = (<9, (j)) . The port k is extended within a 

differential solid angle, dQ.(0,<fi) , at around k = (0,(j)) . We have especially assigned a 
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symbol '0' to the port at a normal direction (0 = 0). The CMT master equation describing 
how a is evolved in time is given by [17] 



da 
dt 



= ico a- 



r i ^ 



V T tot j 



a + x 0 S +0 . 



(l) 



Here co c and 1/ r m (i.e., Q tot = CO c T tot 12) are the resonant frequency and the total decay 
rate of a , respectively. The total electromagnetic energy in the cavity is given by I a I 2 . The 
coupling from the plane wave is described by the last two terms, tc 0 and S +0 . For the 
planewave with an area of A , incident energy per unit time (power) to the 'port 0' is 
(l^+o t ) = = (1 / ^ c£ o n b g l^o T ^ ' wnere s q an d n bg are the permittivity of vacuum and the 
refractive index of the background medium, respectively and ( ) denotes a time average 

I |2 

over one optical cycle. Only a small fraction of the power, LSI , will be coupled to the 




Normally incident 
planewave 



■ |S. 0 I 2 



Photonic-crystal cavity 

Photonic-crystal slab 





Fig. 1. All-dielectric photonic-crystal (PhC) antenna, (a) The 2-D PhC cavity is placed at a 
distance of d from the reflector. The linearly-polarized planewave is being illuminated from the 
top. (b) A coupled-mode theory diagram describes interaction between the incident planewave 
and the PhC cavity. The port k is defined for an imaginary two-way waveguide in the direction 
of k = {6,(j)) with an angular extent defined by a differential solid angle, d€l(6,(j)) . a 
represents the energy amplitude of the cavity. S +k and S. k are the power amplitudes propagating 
in the port k. The PhC antenna is excited only through the 'port 0' with carried power of I S +0 \ 2 . 



resonant mode, which is described by using the coupling constant, k q . According to the 
CMT, k q can be expressed in terms of 1 / t q , the partial decay rate of a only through the 
port 0. Note that although k q is an infinitesimally small value, S +0 and A are infinitely large 

for an ideal planewave. Hence their product form of \k q \ 2 A becomes a certain finite value 
(see Appendix A). 

I |2 

Under continuous excitation, \a\ rapidly grows within the time proportional to ~ Q tot I co c , 
then reaches a certain steady-state value. This saturation energy is obtained from Eq. (1), 
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(\a\ 2 ) I = t 2 \fc\ 2 (\S n | 2 \ = — csn t 2 \e\ 2 • \fc\ 2 A 

\ /rlr^oo tot\ o\ \\ +o\ f T 0 bg tot\ ol I o I 



(2) 



Note that the above equation contains the electric-field intensity of the incident wave, \eJ . 

Another equivalent expression to the above equation can be obtained from the cavity quantum 
electrodynamics (cQED) definition of the mode volume (see Appendix B), 



(H">, 



: -V sn \E. I 2 , (3) 

2 ffl 0 c I k I 7 x ' 



which now includes the maximum electric-field intensity in the cavity. To facilitate an explicit 
expression for |atJ 2 A , we introduce a new term, namely 'far-field coupling factor' [19], 

dP 

(0,4) 

= — — , (4) 

p ot 

which is the ratio between the power radiated per unit solid angle in the direction of 
k = (0,(j)) and the total radiated power, P m . One can easily show that this dimensionless 

quantity varies from 0 to 1 and has the property such that ^ r/(O,0)dQ. = 1 . Furthermore, it 

can also be shown that, for a normally incident planewave, the product term, |/rJ 2 A, is 
directly proportional to rj(6 = 0) (see Appendix A) such that 



kl A = — 

T 



•1/(0 = 0) -cos 2 x- (5) 



Combining Eqs. (2), (3), and (5), finally, we obtain the electric-field intensity 
enhancement factor, 



— 12 



\E I 2 X Q 

®t=T^ = 7 - — -?7(0 = O)-cos 2 z. (6) 

\E\ ti n-n V 

I 0 I bg c m 

This result shows that ® 7 is proportional to both Q tot /V m and rj . Therefore, a PhC cavity 

with the larger Purcell factor [20] is more advantageous for the higher ® 7 . However, for a 
resonant mode that has null vertical emission, rj(6 = 0) = 0 , or an incident planewave with 
X = ti 1 2 , there is no enhancement ( ® 7 = 0 ) regardless of how large Purcell factor is used. It 
may be noted that Q m and rj(0 = 0) are not independent each other because rj(0 = 0) is an 
implicit function of Q m (since Q = co\a\ 2 / P ot ). We also note that most high-g PhC cavities 
(Q m > 100,000) have either null vertical emission or very poor directionality [13]. Therefore, 
we expect a certain trade-off relationship between Q m and r/(& = 0) . In practical 

applications, the maximum Q m may be limited by the bandwidth ( Aa> ) requirement. For 
example, for nonlinear optics applications, where efficient light generation at frequencies 
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different from that of a pumping beam is required, the bandwidth may be chosen to cover the 
wide enough frequency range of interests [6-9]. 
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Fig. 2. Effects of the nano-slot on local-field enhancement. The photonic-crystal (PhC) cavity 

is completely immersed in water ( n = 1.33), where the bottom reflector is excluded, (a) The 

left panel shows a schematic illustration of the nano-slot positioned at the center of the PhC 
cavity. The nano-slot has a rectangular geometry whose width is fixed to 0.05 a while its 

length, L, is to be varied. The right panel shows electric -field intensity (\e\ ) distribution 
around the nano-slot. (b) Quality factor, Q tot , mode volume, V m , and the far-field coupling 
factor, t]{6 = 0) , calculated as a function of L. (c) The electric -field intensity enhancement 

factors, O ? , obtained from the coupled-mode theory (CMT) are compared with those obtained 

from the rigorous finite-difference time-domain (FDTD) simulation. The maximum <£> of 
-750 (-783) is obtained from CMT (FDTD) at L = 0.8 a. 

3. Design of a photonic-crystal optical antenna 

3.1 Effects of a nano-slot 

As a practical example, we have designed a compact, stand-alone, label-free Raman sensor 

[7,21] operating at a wavelength of -1 um. For this purpose, however, Q m should usually be 

limited by -100 in consideration of typical magnitudes of Stokes energy shifts. Designing 
such a low- Q PhC cavity could be challenging as it should have both moderately large 

rj{0 = 0) and small V m . 

Figure 2(a) shows our design, where the conventional single-defect PhC cavity (also 
known as HI) has been modified (see Appendix C). The nano-slot structure introduced in the 
centre of the cavity creates highly-concentrated electromagnetic-fields inside the gap so that it 
can be used to optimize the overlap with an analyte molecule [22,23]. Hereafter, the entire 
PhC cavity is assumed to be immersed in a medium having a refractive index of 1.33 (that of 
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water) for real-time sensing applications. The three key parameters, Q tot , V m , and rj , (see Fig. 

2(b)) can be obtained by using the finite-difference time-domain method (FDTD) [24]. Note 
that far-field emission from our PhC cavity mode is x-polarized at 0 = 0 hence % = 0 in Eq. 

(6). 

First, we consider a PhC cavity without having a bottom reflector. The length of the nano- 
slot, L, is varied while the width of it has been fixed to 0.05 a, which corresponds to -20 nm 

for the operational wavelength (see Appendix D). When L is chosen to be 0.8 a, V m can be 

minimized to be -0.004 A 3 , which shows about 2.2 fold reduction from that of the cavity 
without the nano-slot (L = 0). We note that, as the nano-slot width decreases further, the 
reduction factor for V m can be brought up to -7 times [22] (see Appendix E). Interestingly, the 

condition minimizing V m maximizes Q m so that Q tot I V m can be maximized at the same time. 
Furthermore, rj(0 = 0) in the range of 0.15 -0.18 is already larger than that from the 
conventional dipole- antenna radiation, 3/ 8/r = 0.1 19 (see Appendix F), indicating that the 
efficient inhibition of radiation in the directions parallel to the PhC slab indeed improves the 

free-space coupling. Inserting all the above values, Q m , V m , and rj , into Eq. (6), we get 

considerably large ® 7 of -783, which is already comparable to those obtained from the 

metallic SERS active substrates [25,26]. 

To validate our analytical expression of Eq. (6), we compare our CMT predictions (Fig. 
2(c)) to the rigorous FDTD simulation results that incorporate realistic incoming planewaves 
and subsequent time-evolution of the stored field energy (see Appendix G). As shown in Fig. 
2(c), both results agree well with each other indicating that our CMT formalism is valid even 
for a resonance with Q < 100. It should be noted that, however, it may not be applicable to the 
conventional metallic antennas, because they usually support too low Q (typically, <10) [1- 
4,25,26] to justify the weak-coupling assumption of CMT [18]. 

3.2 Effects of a bottom reflector 

We now consider the effect of a bottom reflector on far-field radiation patterns. In the absence 
of the bottom reflector, radiation losses occur symmetrically with respect to the PhC slab as 
shown in Fig. 3(c). If the downward-emitting components could be suppressed, one would 

double rj(6 = 0) as well as the final intensity enhancement, ® 7 (see Eq. (4)). In fact, the gap 

distance, d, is a key parameter to control such enhancement through interferences between the 
originally upward- and downward-emitting waves, both of which can be met at 6 = 0 by the 

reflector (see Fig. 4a) [16]. As a result, both far-field radiation (or rj(6 = 0) ) and Q m can be 

modified. For example, as shown in Fig. 3(d), if d equals to the effective 1.0 X distance 

( d = A I n bg = A 1 1.33 ), fairly-good unidirectional beaming can be obtained [16]. In this 

regard, our design is reminiscent of the Yagi-Uda antenna [5] in the radio-frequency domain. 
What is interesting here is that one could achieve very high r/(@ = 0) -0.67 at 

d = 0.95/1 / n b (see Fig. 4(b)), which is about 4 fold enhancement in comparison with that of 
the PhC cavity without a bottom reflector. Strong modulations in Q m , V m , and rj as a 
function of d with a periodicity of - 0.5 x (A I n bg ) are clearly visible, indicating that the far- 
field interference is associated with the gap size [16]. 
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Fig. 3. Far-field radiation patterns. Far-field patterns are calculated by using the finite- 
difference time-domain method (FDTD) and the near- to far-field transformation formulae 
presented in ref. 11. The two far-field patterns (x,y) in the upper row are for the upper 

hemispherical surface (0 < 6 < 90 ), where the mapping between (x,y) and ( 6 , $ ) is given 
by x = 6 cos (ft and x = Osincj) . In the two polar plots in the lower row, the dotted line (solid 
line) shows the radiation pattern in the x = 0 (y = 0) plane, (a, c) The photonic -crystal cavity 
without the bottom reflector, where the length of the nano-slot is 0.8 a. (b, d) The same 
photonic-cavity cavity as before but with the bottom reflector, where the gap size is 

d = X I n, . 

bg 

Similar comparison between our CMT results and rigorous FDTD results has been 
performed for the PhC antenna with the bottom reflector. For d < 0.9 A I n bg , the both results 

seem to agree well with each other. However, as d increases further, the CMT prediction 
begins to fall below the exact FDTD result. The CMT model predicts that the maximum of 

0 7 will occur at d = 0.9A/n bg with ® 7 -3,437. Clearly, the PhC optical antenna 

outperforms or at least shows similar performances compared with metallic nano-structures 
[1-6,8,9,25,26]. However, the exact FDTD simulation predicts higher enhancement of -4,500 

at a slightly different gap size of d = 1 .0/1 / n b . Such discrepancy between the two 

predictions might be due to the fact that the gap region may act as a Fabry -Perot like cavity 
enclosed by the top PhC slab mirror [16]. This broad resonance may interfere with the narrow 
resonance of the PhC cavity to produce a Fano-like spectral feature [17]. The result shown in 
Fig. 4(c) suggests that the more rigorous CMT model incorporating all those interactions 

should be developed for d > ~ A n b . 
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Fig. 4. Effects of the bottom reflector on local-field enhancement. The length of the nano-slot is 

fixed to be 0.8 a, at which O has been optimized in the absence of the bottom reflector, (a) 

Electric-field intensity distribution in the x-z plane. The downward-propagating waves can be 
redirected to the top by the reflector and they can interfere with the originally upward- 
propagating waves. The interference condition can be determined by the gap size, d. (b) 
Quality factor, Q tot , mode volume, V m , and the far-field coupling factor, rj{6 = 0) , calculated 

as a function of d. (c) The electric -field intensity enhancement factors, <I> ; , obtained from the 
coupled-mode theory (CMT) are compared with those obtained from the rigorous finite- 
difference time-domain (FDTD) simulation. <D of -3,437 (-4,500) is obtained from CMT 

(FDTD) at d = 0.9 A I n bg ( d = 1.0 A I n bg ). 

4. Concluding remarks 

The proposed design still leaves room for further improvement. For example, by narrowing 
the width of the nano-slot, additional factor of 2-3 times <D 7 is achievable. By employing 

high-g (Q tot > 10,000) modified dipole modes, ® 7 > 10 5 is expected from our CMT model 

(see Appendix H). For enhancing nonlinear signal generation based on such high- Q antennas, 
one may utilize the multitude of resonances; one for the pumping wavelength and others for 
the nonlinear-converted wavelengths [27]. Note that such dramatic results do not rely on 
technologically challenging geometries, such as extremely sharp tips and/or small gaps, which 
are common features in many metallic antenna designs [1-3,6,19]. Furthermore, one may 
develop a miniaturized antenna for free- space communication. The similar design can also be 
applied to high-efficiency photovoltaic devices for driving nano-electronics [28]. The 
proposed scheme is fully compatible with the matured silicon nanofabrication technique and 
could be used in a wide range of applications where strong light- and-matter interactions are 
essential. 
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Appendix A - Calculation of |at 0 | A. 

According to the known result from Fraunhofer diffraction theory, the finite extent (area A ) 
of an incoming planewave inevitably involves the angular spread, AQ , such that 

(A/n ) 2 

AQ = ^- (A-l) 

A 

where X is the wavelength in vacuum. From the well-known result of CMT [17], we obtain 
the relation between the coupling constant rc Q and the partial decay rate 1 / r 0 such that 

|atJ 2 = (2/ ?" 0 )cos 2 x ' wnere X * s me an gl e between the incident wave's polarization and the 

outgoing radiation's polarization. The ratio between the power radiated into the solid angle, 
AQ , and the total radiated power is given by 

AP(0*0) _1/t q 



P 1/r 

tot to 

Combining all the above information, 



I I 2 . 2 

kl A = - 



1 2 2 

cos x = — 

AQ r 



f „ A 



\ n H) 



(A-2) 



AP(fl«0) 1 2 

cos x (A-3) 

AQ P 



In the limit of infinitely large area, as both AP and AQ become infinitesimally small, the 
above formula equals to Eq. (5). 

Appendix B - Definition of mode volume V m . 

The cQED mode volume is defined as the ratio of the total electromagnetic energy stored in 
the cavity and the peak electromagnetic energy density, such that 

jjj-s 0 s(?)\E(rJ)\ 2 +- Mo \H(rj)\ 2 d^ 
V = ~, ~ r, (B-1) 



max e 0 e(f) |£(F, t)f + ^ ju Q \li(r 9 1)\ J 



where we have assumed that the cavity is consisting of non-magnetic materials (ju(r) = 1 ). 
Taking time averages on both the denominator and the numerator and using the fact, for 

harmonically oscillating fields, (^(1/ 2)s 0 s\E\ 2 ^j = (^(1/ 2)jU 0 \h\ 2 ^ , one can rewrite the 
above expression as 

(kl 2 ) (kl 2 ) 

V =-. J—* '-i 7T^ = ^ — (B-2) 



(max{f 0 f(?)|£(F,f)| 2 }) r 



1 P 

^ 0 c I loc I 



where n and E 1 are the refractive index and the electric-field at the location where the term 

c loc 

braced is maximized. 
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Appendix C - Design of the dipole-mode cavity with Q < 100 



A single-defect PhC cavity is formed in a 2-D triangular-lattice PhC slab. The slab thickness, 
T, and the refractive index of the slab, n s i ab , are chosen to be 0.6 a and 3.5 (that of Si at near- 
infrared wavelengths), respectively, where 'a' denotes the lattice constant of the PhC. For 
potential real-time bio-sensing applications, the refractive index of a background material is 
assumed as water, so n bg = 1.33. For a triangular air-hole type PhC slab, the radius of air-holes 
is usually chosen to be R ~ 0.35 a, which allows mechanical robustness as well as wide- 
enough PBG in the in-plane directions. Throughout the paper, we fix R to be 0.36 a, which 
gives the in-plane PBG from co n = 0.27 ~ 0.35, where oo n is the normalized frequency, a/A,. 

However, without introducing any structural modification, the cavity supports doubly- 
degenerate dipole-modes (the two dipole-modes are at the exactly same frequency) due to the 
perfect six-fold symmetry. In order to break the inherent degeneracy and to lower the Q a bit 
more, we have changed positions of the four-holes as a function of their radii (R m ) such that 
(see Fig. 5(a)) 



Modified holel: 



Modified hole2: 



Modified hole3: 



a — » 



v 2 2 y 



4i 



V 2 2 J 



a — » 



a -> 
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Modified hole4: 



V 
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a — » 



J 



V 
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2 2 



We have investigated how Q and oo n are evolved as we vary R m (Fig. 5(b)). Maximizing 
the strength of in-plane confinement is desirable for higher rj{6 = 0) . Therefore, we choose 

R m to be 0.36 a, which results in Q of -45. For the deformed dipole mode (R m = 0.36 a), 
electric-field intensity distribution (|^| 2 ) as well as its electric-field lines are shown in Fig. 

5(c). We can see that it has the highest \e\ 2 around the centre of the cavity, where electric- 
fields are mostly parallel along the x direction. This is the main reason why the dipole-mode 
has strong vertical radiation with the well-defined linear polarization (x-pol) (see Fig. 5(d)). 
This fact also provides us an additional hint for further optimization of the mode volume, 
which will be explained in Appendix E. It is a known fact that the other degenerate- split 
dipole pair is orthogonally polarized (y-pol) to this dipole-mode [11]. 
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Fig. 5. (a) Structure of the \ow-Q dipole mode cavity, (b) Normalized frequency and Quality 
factor as a function of the modified hole radius, R m . (c) Vector plot of the electric -fields, (d) 
Far-field radiation pattern from the mode shown in Fig. 5(c). 



Appendix D - Tuning of the PhC cavity by varying the length of a nano-slot 

Figure 6 shows evolution of the normalized frequency (oo n ) and the resonant wavelength (in 
case of a = 350 nm) as a function of the slot length (L). In this paper, the width of the nano- 
slot is fixed to 0.05 a, which corresponds to 17.5 nm. When L is varied between 0.6 a and 1.0 
a, the resonant wavelength is quite stable at ~ 1,095 nm. Slight tuning the lattice constant by 
about 2% will bring down the resonant wavelength to 1,064 nm, which is suited for Nd:YAG 
laser pumping. 



= 350nm assumed 




0.315 

D.O 0.4 0.8 1.2 1.6 

Slot Length (a) 

Fig. 6. Resonant wavelength (nm) and normalized frequency (a/A,) as a function of the slot 
length (a), where a = 350 nm was assumed. 

Appendix E - Introduction of a nano-slot into a PhC cavity 

As shown in Fig. 2(a), we introduce the nano-slot in the region where electric -fields are 
aligned along the x direction. If the width of it is sufficiently small not to perturb the original 
field distribution too much, then we can safely assume that the electric -fields are still aligned 
in the same direction. However, the amplitude of the electric-field within the slot can be 
changed noticeably due to the boundary condition for the normal components of the electric- 
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fields. The similar idea was proposed by J. T. Robinson, et al. and applied to 1-D nanobeam 
cavities [22]. Here, we briefly discuss how the nano-slot can reduce the mode volume. 

First let us recall the boundary condition for normal components of the electric-field on 
either side of the boundary surface (at x = ±A/ 2 ), 

n lh E (x = A/2) = n h E (x = A/2). (E-l) 

slab x v 7 bg x v 7 v/ 

Here E x and E x are the electric-fields before and after introducing the nano-slot at the centre 

of the cavity, respectively. We also assume that the width of the slot is so small ( A — » 0) that 
it could not change the total electromagnetic energy contained in the cavity. From the 
definition of the cQED mode volume, before introducing the nano-slot, we have 

<H 2 > 

y = - i — . (E-2) 

- f o«Lk(* = °)f 

After introducing the nano-slot, however, 



v. ¥\ 



( n ^ 



1 2 I . \: 
e o n J E S x = °)\ 



y . (e-3) 



2 ° b8 



Therefore, for an ideal nano-slot ( A « 0 ), the maximum reduction in the mode volume could 
be up to (n lab I n bg ) 2 = (3.5 / 1.33) 2 « 6.9 . 

Appendix F - Examples of the far-field coupling factor, t\ 

For isotropic radiation pattern over the entire An steradian (ideal spherical source) 

r dP dP r dP 
Because dP I dCl is a constant, P = dQ, = \dQ = 47r . Therefore, we get 

tot J dn dn : da 

7](0, (j)) = rj(0 = 0) = 1 / 4tt = 0.0796 

For classical dipole antenna radiation 

dP I dCl is given by dP I dfl = (3P ot 1 8;r)sin 2 (;r 12-6), where the direction of the dipole 
oscillation is assumed to be in the x-y plane. Therefore, rj(6 = 0) = 3 / %n = 0. 1 19 . 

Appendix G - Calculation of 0 7 using FDTD 

In order to estimate the electric-field intensity enhancement (O y ) through Eq. (6), the three 
important parameters, Q,V m , and rj(6 = 0) , should be obtained by using FDTD. 

Cavity Q 

E x -field component at the center of the cavity has been detected as a function of time. Then 
the harmonic inversion algorithm (harminv) [29] is used to obtain Q. 
Mode volume V m 

For sufficiently high- Q modes, both of the following definitions will produce the almost same 
values. 
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\\\s 0 s(?)\E(rjj\ 2 d 3 r 
V = AAJ ~n — —pr , (G-l) 



max|^(r)|i(r,^)| 2 } 



V (2) = ^— ^ ^ (G-2) 



(max{^(F)|i(F,0| 2 }) r 



Here, in the first definition, t m is the time when harmonically varying electric-fields are 
maximized. In the second definition, time averages are performed over one optical period. 
However, for low-g modes, differences between and V w (2) become more pronounced 
because of the fast decay of the fields within one optical cycle. We have used the second 
expression, V^ 2) , for our definition of the mode volume (see also Appendix B). 

Far-field coupling factor, r/(0 = 0) 

As depicted in Fig. 7(a), we set up a small area (S 0 = (a/40) x(a/40) ) which locates 
sufficiently far away from the cavity (L z > 3X). Then, the power radiated through the whole 
surface (S tot ) and the power radiated through the small area (So) are obtained during FDTD 
time stepping. Using the definition of rj(6 = 0) (see Eq. (4)), 



7/(0 = 0) : 



dP 




d£l 


0=0 


P 



f Power radiated through S ^ 

& 0 



S IE 

tot - 0 z 

Rigorous FDTD simulation for O y 



1 



Power radiated through S 



(G-3) 



As depicted in Fig. 7(b), a continuous-wave (CW) planewave source (\E Q \ ) which is 
polarized along the x direction is generated and propagated onto the PhC cavity. We detect 
\E t oc (t)\ 2 as a function of time at the centre of the slot until it reaches a certain steady-state 

value. Then, we compute the ratio between (^E hc (t °°)| 2 ^ and (|^ 0 | 2 ) t0 S et me 

enhancement factor, O 7 . 
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<WH> 




Fig. 7. (a) FDTD simulations used to obtain the far-field coupling factor, rj. (b) The rigorous 
FDTD simulation used to directly obtain O . 

Appendix H - Extremely large local-field enhancement > 10 s from the modified dipole 
mode 

We have designed a high- Q dipole mode by modifying the six nearest neighbor holes as 
depicted in Fig. 8(a). The slab thickness, T, and the refractive index of the slab, n s i ab , are the 
same as before. The entire device is assumed to be completely immersed in water (n bg = 1.33). 
The radius of background air-holes is chosen to be r = 0.35 a and the six nearest holes are 
slightly pushed away from the cavity center after reducing their radii to r m = 0.28 a. To 
estimate resultant local electric-field intensity enhancement, we have used Eq. (6). Three 

important parameters, Q , V m , and rj(0 = 0) , along with the resultant field enhancement 

values are presented in Table 1 for several representative cases with and without a bottom 
mirror. 



Table 1. Electric-field intensity enhancement from the high-g dipole mode. 





Without a 
mirror 


Effective gap 
0.95 X 


Effective gap 
1.00 X 


Effective gap 
1.05 X 


Q t o t 


5812 


5464 


5261 


5256 


V m (X 3 ) 


0.00315 


0.00316 


0.00315 


0.00315 


i/(e = o) 


0.1956 


0.2267 


0.2989 


0.2850 


\E loc \ 2 AE 0 \ 2 


97,682 


106,246 


134,948 


128,632 



It is noteworthy that rj(6 = 0) has been compromised with the high-g, to a relatively low 

value of -0.3. However, owing to their high Q/V values, one can achieve extremely high 
electric-field intensity enhancements over 100,000. Highly localized electric-field intensity 
distributions within the nano-slot are shown in Fig. 8(b) and Fig. 8(c). 
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Fig. 8. (a) Structure of the high-<2 PhC dipole mode cavity. (b,c) Electric -field intensity (IEI 2 ) 
distributions (b) in the x-y plane and (c) in the x-z plane. 
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